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Abstract. We show that the dynamics resulting from preparing a one- 
dimensional quantum system in the ground state of two decoupled parts, 
then joined together and left to evolve unitarily with a translational invariant 
Hamiltonian (a local quench), can be described by means of quantum field theory. 
In the case when the corresponding theory is conformal, we study the evolution of 
the entanglement entropy for different bi-partitions of the line. We also consider 
the behavior of one- and two-point correlation functions. All our findings may 
be explained in terms of a picture, that we believe to be valid more generally, 
whereby quasiparticles emitted from the joining point at the initial time propagate 
semiclassically through the system. 
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1. Introduction. 

The understanding of the degree of entanglement in extended quantum systems has 
prompted an enormous effort in a field at the border of condensed matter physics, 
quantum information theory, and quantum field theory. It would take too long to 
mention all the topics and problems addressed so far: we refer the reader to a recent 
thorough review and extensive bibliography [1] . 

Several measures have been introduced to quantify the entanglement in an 
extended system [JJ. Among them we consider here only the entanglement entropy, 
which is defined as follows. Suppose a whole system is in a pure quantum state 
with density matrix p = |\&)($|, and an observer A measures only a subset A of a 
complete set of commuting observables, while another observer B may measure the 
remainder. A's reduced density matrix is pa = Trap. The entanglement entropy is 
just the von Neumann entropy 

S A = -Tta Pa log pa (1) 

associated with this reduced density matrix. For an unentangled product state, 
Sa = 0. Conversely, Sa should be a maximum for a maximally entangled state. 
One of the most striking features of the entanglement entropy is its universal behavior 
displayed at and close to a quantum phase transition. In fact, close to a quantum 
critical point, where the correlation length £ is much larger than the lattice spacing 
a, the long-distance behavior of the correlations in the ground state of a quantum 
spin chain are effectively described by a 1+1 dimensional quantum field theory. At 
the critical point, where £ diverges, the field theory, if Lorentz invariant, is also a 
conformal field theory (CFT). In the latter case, it has been shown that, if A is an 
interval of length I in an infinite chain, Sa — (c/3) log£ [2j [3J H] , where c is the central 
charge of the corresponding CFT. When the correlation length £ is large but finite 
(more precisely £ -C £), it has been shown that, increasing £, Sa saturates [3] to 

s A ~(c/3)io g £ mm. 

Recently the interest in the properties of entanglement has been extended to 
understanding its dynamical behavior. The natural question is how the entanglement 
propagates through the system when it is prepared in a state that is not an eigenstate 
and then is left to evolve in the absence of any dissipation. The most common situation 
studied so far concerns a sudden quench of some coupling of the model Hamiltonian 
O [7] [H [10] . The main motivation for the large number of studies of this dynamics 
is that it may actually be realized and measured in ultra-cold atomic systems 
For a complete list of references on the subject we refer to our earlier paper [12] . 

Here we consider a different situation, known in the literature as local quench 
[13] , We will concentrate on the case when a physical one-dimensional system (e.g. a 
spin chain) is physically cut into two parts that are in their own ground state. We 
then join the two halves at a given time t = and study the subsequent evolution 
according to a translational invariant hamiltonian. Here, the entanglement entropy is 
the most natural quantity to be studied because the two halves are clearly unentangled 
for t < whereas the initial energy differs from that of the ground state only by a 
finite amount. This differs from the case of a global quench, when the energy of 
the initial state above the ground state is always extensive, and correspondingly the 
entanglement can be extensively larger. We will only consider gapless models that are 
described asymptotically by a boundary CFT, since in this case we can use powerful 
analytic results [HI [15] . 
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The paper is organized as follows. In Sec. we briefly recall the CFT approach 
to the entanglement entropy. We outline the setup for the local quench which is then 
applied to different situations in the following section [31 Then we consider the time 
evolution of correlation functions in Sec. [U In Sec. [5] we describe how this approach 
can be generalized to the case with two joining points and we derive an heuristic 
solution. Finally we conclude the paper with a discussion of open problems and 
possible generalizations in Sec. El 



2. The CFT approach to entanglement entropy and a local quench 

2.1. Entanglement entropy and CFT 

The entanglement entropy Sa defined by Eq. ([TJ) can be studied in a general quantum 
field theory through the replica trick [3] 



(2) 

=i 



This is particular useful in a CFT because, when A consists of disjoint intervals with N 
boundary points with B, Trp^ transforms under a general conformal transformation 
as the N point function of a primary field with conformal dimension [4] 



c 1 



12 \ n n 



(3) 



where c is the central charge of the underlying CFT. 

In particular this implies that the entanglement entropy of a slit of length £ in an 
infinite system is given by 



A 2x " c, £ 



Trp J X = c n - S A = -log-+4, (4) 

\a ) da 

where a is an UV cutoff (e.g. in a spin chain is the lattice spacing). The constants c ra 
are also non-universal and so is the derivative for n = 1, c[, entering in Sa- However, 
in any particular model, a and c[ can be fixed unambiguously by specific requirements 
(see e.g. [5]). 

Another important result we will use in the following is the entanglement of the 
slit A = [0,£] in a semi- infinite line with specific boundary conditions at r = 0, that 
is given by [4] 

Trpl = c n (-) ^ => S A =%g-+c' 1 . (5) 
\ a J b a 

In analogy to the c„, the constants c n are not universal and depend on the boundary 
condition at r = 0. However comparing the finite temperature result of Ref. [3] 
with the standard thermodynamic entropy for systems with boundaries [16] we have 

01 El Eg 

c' 1 -c' 1 /2 = logg, (6) 

where logg is the boundary entropy first introduced by Affleck and Ludwig [16] . We 
recall that g is universal and and depends only on the boundary CFT. 
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2.2. CFT for global quenches 

CFT is also able to predict the time-dependence of the entanglement entropy after a 
global quench [B]. We briefly summarize here these results, because some of them will 
be useful in understanding the local case. The density matrix has the path integral 
representation [6] 



W(r")\p(tW(r')) = Z-\r(r'')\e- UH -' H \Mr)){Mr)\e +UH -' H \^\r')) , (7) 



where we included damping factors e~ in such a way as to make the path integral 
absolutely convergent. Each of the factors may be represented by an analytically 
continued path integral in imaginary time: the first one over fields ^(r, r) which take 
the boundary values 4>o{r) on t = —e — it and 4>"{ r ) 011 r = 0, and the second with 
ip(r,r) taking the values tp'(r) on r = and V'oM on r = e — it. When H is critical 
and the field theory is a CFT, under the renormalization group any translationally 
invariant boundary condition is supposed to flow into a boundary fixed point, satisfying 
conformal boundary conditions. Thus we may assume that the state \ipo) corresponds 
to such boundary conditions on sufficiently long length scales. Thus, for real r, the 
strip geometry described above may be obtained from the upper half-plane (where the 
entanglement entropy is known from [4]) by the conformal mapping w = (2e/7r) logz. 
In this way one easily get the result for imaginary time. The real time evolution is 
obtained by continuing r = e + it. The final result is [5] 



From this we note that e enters in the calculation in an essential way. A careful analysis 
[19l \12\ shows that in fact it corresponds to the correlation length in the initial state 
and so it is not only a useful tool, but a physical important parameter. 

The fact that 5a (i) increases linearly until it saturates at t* — 1/2 has a simple 
interpretation in terms of quasiparticles excitations emitted from the initial state 
at t = and freely propagating with velocity v = 1. This phenomenon is very 
general and holds even for non-critical systems, and it has been confirmed by exactly 
integrable dynamics and numerics (see e.g. 0E1IT3]) However, since in lattice models 
there are particles moving slower than v, after t* the entropy does not saturate 
abruptly, but is a slowly increasing function of the time. The same picture is valid 
for the correlation functions. Firstly incoherent quasi-particles arriving a given point 
from well-separated sources cause relaxation of (most) local observables towards their 
ground-state expectation values. Secondly, entangled quasiparticles (emitted from an 
initially correlated region) arriving at the same time t induce correlations between 
local observables. In the case where they travel at a unique speed v, therefore, there is 
a sharp "horizon" effect: the connected correlations do not change significantly from 
their initial values until time t ~ |r|/2u. After this they rapidly saturate to time- 
independent values. For large separations, these decay exponentially differently from 
the power law dependence in the ground state. Also for correlation functions, this 
picture has been shown to be valid in several models (see e.g. flUl [T2l [2131 |2"TI l2"2l |2"5] ) . 

It is clear that a very similar interpretation should apply also for local quenches, 
but with the fundamental difference that now the quasi-particle excitations are emitted 
only from the point where the quench happened and not from everywhere. 




(8) 
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Figure 1. Space-time 
region for the density 
matrix (left) mapped to 
the half-plane (right) by 
means of Eq. {5}. z\ = 
ir and 22 = «T+f during 
the computation and in 
the end r — > it. 



,2.3. CP T setup for local quenches 

Suppose we physically cut a spin chain at the boundaries between two subsystems A 
and B, and prepare a state where the individual pieces are in their respective ground 
states. In this state the two subsystems are completely unentangled, and its energy 
differs from that of the ground state by only a finite amount. Let us join up the pieces 
at time —t and watch the system evolve up to t = 0. 

The procedure for the global quenches does not apply because the initial state is 
not translational invariant and will not flow under the renormalization group toward a 
conformally invariant boundary state. One may try to handle the problem introducing 
proper boundary condition changing operators (see e.g. |24j ) . However we prefer 
to take a different approach. In fact, we can simply represent the corresponding 
density matrix in terms of path integral on a modified world-sheet. The physical cut 
corresponds to having a slit parallel to the (imaginary) time axis, starting from — oo up 
to 7~i = — 6 — it (the time when the two pieces have been joined), and analogously the 
other term of the density matrix, like in Eq. ([7]), gives a slit from r 2 = e — it to +oo. 
Again we introduced the regularization factor e, that we will interpret a posteriori. 

For computational simplicity we will consider the translated geometry with two 
cuts starting at ±ie and operator inserted at imaginary time r. This should be 
considered real during the course of all the computation, and only at the end can 
be analytically continued to it. This plane with the two slits is pictorially represented 
on the left of Fig. [T] where ir corresponds to z\. As shown in the same figure, the 
z-plane is mapped into the half-plane Re w > by means of the conformal mapping 

z If w — 1 

w — — hi/ - +1 with inverse z — e . (9) 

e V Ve/ 2w 

On the two slits in the z plane (and so on the imaginary axis in the w one) conformal 
boundary conditions compatible with the initial state must be imposed. For example, 
in the most natural situation when the boundary "spins" are left free, we require free 
boundary conditions. Oppositely, when the boundary spins are supposed to stay in a 
particular state we require fixed boundary conditions. 

Note that this setting is valid in any dimension when the system is prepared in 
two spatially divided halves. In this case one can try to tackle the problem with the 
methods of boundary critical phenomena |25j . but this will be extremely difficult, if 
not impossible, since conformal invariance is far less powerful. 



Entanglement and correlation functions following a local quench: a CFT approach 6 



D 
ID 
III) 
IV) 



A 



f 



Defect 



_B_ 

B 



B 



B 



Figure 2. The four different bipartitions of the line we consider here. 



3. Entanglement entropy 

In this section we consider the time evolution of the entanglement entropy after the 
local quench. We consider two half-chains joined together at the point m — at a 
given time. The various subcases we consider correspond to different spatial partitions 
of the system among which we calculate the entanglement. We consider the four 
different situations depicted in Fig. [2] 

We start with the more natural division, considering the entanglement entropy 
between the two parts in which the system was divided before the quench (case I). As 
case II we consider the entanglement of the region r > £ with r < I. This will allow 
us to identify simple horizon effects. To highlight the interaction of the two previous 
effects we then consider the entanglement of the part B = [0, £] with the rest (case 
III). Finally we consider as IV the most general slit B = [ii,^]- Clearly cases I and 
III are particular choices of II and IV respectively and can be derived as simple limits. 
However we prefer to present the results in such order because the physical effects and 
their interpretation will emerge in a more natural way 



3.1. Case I: Entanglement of the two halves 

This is the case when B is the positive real axis and A is the negative real axis. Tr p\ 
transforms like a one-point function that in the w plane is [2Reu;i] _x " . Thus in the z 
plane at the point Z\ = (0, it) we have 



that using 



becomes 




($n) = C„ 



ae/2 



^ € — T 

Continuing this result to real time r — > it we obtain 
(1 , _ / ae/2 

(*»(*)> = c » [jrfT2 

Using finally the replica trick to find the entanglement entropy we have 



t 2 +e 2 



= — log 

6 B ae/2 



(10) 

(11) 
(12) 
(13) 
(14) 
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There are two main pieces of information that we can extract from this result, coming 
from the long and the short time behavior. For f>£we have 

S A (t»e) = flog -+**,, (15) 
6 a 

i.e. the leading long time behavior is only determined by the central charge of the 
theory in analogy with the ground state value for a slit. This could result in a quite 
powerful tool to extract the central charge in time-dependent numerical simulations. 
This behavior resembles the slow increasing of Sa observed numerically in Ref. [26] 
(unfortunately these data were never fitted, to understand the quantitative predictive 
power of our formula). The constant k is given by ko = c\ + (c/6) log(2a/e). 

The behavior for short time allows instead to fix the regulator e in terms of the 
non- universal constant c^. In fact we have 

S A (t = 0) = Uog-+c' 1 =0 => e =V 6S 'i/ c . (16) 
o a I 

Note that, even if non-universal, for a given lattice model the constant can be 
fixed only by ground state quantities. Consequently Eq. (|14p has no free dynamical 
parameter. 

The parameter e depends in a specific manner on the boundary contribution to 
the entanglement c^. This is completely different from what happens for a global 
quench where the equivalent regulator is connected to the correlation length in the 
initial state [HI [12] . Furthermore the order of magnitude of e is fixed by the lattice 
spacing a (i.e. the UV cutoff) and so considering the limit of times and lengths larger 
than e is equivalent to the standard condition to apply the field theory for distances 
larger than a. 



3.2. Case II: Decentered defect 

Let us now consider the entanglement of the region r > I with the rest of the system. 
In this case Tr p n A is equivalent to the one-point function in the plane z at the point 
Z2 = t + ir as in Fig. [1] Under the conformal mapping ([1]) this point goes into 



ew 2 



= i + ir + ^e 2 + (£ + ir) 2 = I + ir + 



pe 



with 



p 2 = ^2 + (2 _ r 2)2 + U 2 T 2 

and 



1 2£t 
— arctan — ■ 

2 Ve 2 + 1 2 



p COS ( 



p sm ( 



dw 
dz 

Thus ($„) is given by Eq. (fT2"|) with w% — > W2, resulting in 



ay/(£ + p cos 0) 2 + (t + p sin t 



2p(t + pcos6) 



(17) 
(18) 

(19) 
(20) 



The real time evolution for t, I 3> e is obtained by analytically transforming r — > it. 
The calculation can seem very cumbersome, but it is greatly simplified by the fact 
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that in this limit we have p 2 — > \£ 2 — t 2 \ and pcosO — > max[£,t], ps'm9 — > imm[£,i\. 
Care must also be taken when the zero-th order in e is vanishing. In the end we have 

i 1 

y/(£ + pcos9) 2 + (T + psm9) 2 J 2|£| 



2p(£ + pcos6) 



2{t 2 - I 2 ) 



t <£, 
t >£, 



(21) 



that leads to 



S A 



c, 21 _ 
- log - + c[ 
6 a 



c ^ 
- log — 

6 B 



t <£. 



t > 



(22) 



with ko the same as in Eq. (|15p . The interpretation of this result is quite direct. 
Indeed at t = the joining procedure produces a quasi-particle excitation at r = 
that propagates freely with the corresponding speed of sound v s that in the CFT 
normalization is v s = 1. This excitation takes a time t = £ to arrive at the border 
between A and £? and only at that time will start modifying their entanglement. The 
following evolution for t 3> £ is the same as in Eq. (Ti"5j) . 

Also the constant value for t < £ deserves a comment: it is exactly the value 
known from CFT for the slit in the half- line Eq. {5} • This is a non-trivial consistency 
check. Note that a finite e smooths the crossover between the two regimes and makes 
the entanglement entropy a continuous function of the time. 



3.3. Case III: The slit with the defect at the border 

Let us consider again the same physical situation as before, but we now calculate the 
entanglement entropy of A = [0, £] and B the remainder. For t < the real negative 
axis is decoupled from the rest and does not contribute to the entanglement entropy 
that is just the one of a slit in half-chain, i.e. the initial entropy is given by Eq. (0. 

The entanglement entropy is obtained from the replica trick considering the 
scaling of a two-point function between the endpoints of the slit. In the z plane 
these two points are z\ = it and zi — £ + ir (we adopt the same notation as before to 
make direct the use of previous formulas). As usual, the plane with cuts is mapped 
into the half-plane Kew > by the transformation ([9]). On the w plane we have [4] 



a \Wi + W2\\W2 + W\\ 



($ n («;i)*_ n (tD2)) = C 2 n 

\ |lOi - w 2 \\w 2 - W\\\W\ + Wl\\Wi 

so the mapping to the original plane z is 



■ U> 2 



dw 



dz 



dw 



dz 



a^\wi + w 2 \\w 2 + Wi\ 



\Wi - w 2 \\w 2 - Wx\\Wl + Wi\\W 2 + w 2 \ 



(23) 



(24) 



e 2 K 
e 2 |wi 



where the various terms are functions of z; through Eqs. (TTTj) . (jTTJ) , (fTS]) , (jT5J) , and 

W2 \ 2 = (£ + pcosO - Ve 2 -t 2 ) 2 + p 2 sin 2 6, (25) 
- W2 1 2 = (£ + p cos 9 + \/<l 2 - t 2 ) 2 + p 2 sin 2 9 . (26) 
Putting everything together we get 
($ n (zi = iT)$- n (z 2 =ir + £)) = 
pcost 



-2\2 



p 2 sin 2 i 



{£ + p cos 6 - Ve 2 



r 2 ) 2 + p 2 sin 2 I 



y/{£ + pcosO) 2 + (r + psinf 
Ap{£ + p cos 9) 



(27) 
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Again this simplifies considering the limit t, £ e as explained after Eq. (I^D|) . We 
finally obtain 



T'7>",H (28) 

Note that for t > £ the calculation is slightly more cumbersome because we should 
take into account the first order in e of some expressions. 

Using finally the replica trick we get for the entanglement entropy 

c , t c , £ c . ,l — t 
ln- + -ln- + -In 4 j— + K t<£, 

Sa={ 3 c \ 6 6 6 £ + t (29) 

% In - + 2c\ t>£. 
S a 

The crossover time t* = £ is again in agreement with the quasi-particles interpretation. 

There are several interesting features of this result. For very short time t <C £ it 
reduces to the £ ~ oo case Eq. (IT5t as it should. The leading term for t > £ is just 
the ground state value for a slit in an infinite line. However the subleading term is 
not the same, in fact we have 

S A (t>£)- S A (eq) = 2ci - ci = log g , (30) 



where we also use the value of e in Eq. (|16|) . This is a signal that for long time 
the system still remembers something of the initial configuration as a boundary term 
that is unable to "dissipate" . Since the extra energy never dissipates under unitary 
evolution, there is no reason for the constant terms to be the same. 

Another interesting feature is the behavior for t < £. This is very similar to the 
form proposed in Ref. |13j to fit the numerical data, i.e. 

S A = C f \og£+ | log(t^) + | log(l - t/£) + k' . (31) 
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Only the term in t + £ was missing in Ref. |13] . However this behaves smoothly for 
< t < £ and its effect can be well approximate by a constant factor in k' . It will be 
interesting to check whether and how the use of this term in £ + 1 changes the quality 
of the fit. Furthermore the fit of Ref. [T3] is also a strong confirmation of our result, 
indeed they found cq ~ 1 + C2, ci ~ 1, and c 2 ~ 1/2 (plot 7 in Ref. [13] with t' = 0) 
that are exactly our predictions for c = 1 (the known central charge of the XX model) . 

Furthermore Eq. (|29p display a large plateau for 0.2 < t/l < 0.8 as shown in Fig. 
[3] (the black line in both the plots) and already noticed in Ref. [13] . This plateau can 
be studied considering the value at the maximum t/£= (V5- l)/2 where we have 

c I 3 

Sa (plateau) = - log - + h = -S A (cq) + k 2 , (32) 
2 a 2 

with the constants k p simply related to a, e, and c\. 

Finally we can study the result for t = when e -C £■ 

c Ae£ c 21 

S A (t = 0) = - log — + 2c[ = - log - + c[, (33) 
o a ba 

where in the last equality we used the value of e given by Eq. ([T6]) . This result is 

exactly the initial value for a slit of length £ in the half-line. Again this is an highly 

non-trivial consistency check of all the method. A final curiosity is that 

2S^(t = 0)-S A (t>0 = |log2, (34) 

is independent of the details of the theory. 



3.4- Case IV: The general slit 

Let us now consider the most general case of a slit A — [£2,(1]- In the z plane we have 
to calculate the two-point function between z\ = £\ + ir and z 2 = £2 + it. Under the 
mapping (p} these correspond to w p = l p + ir + p p e l6p with p p and 9 p given by Eq. 
(|18|) with £ — > £ p and p = 1,2. In imaginary time, after a long, but straightforward 
algebra we have 

($ n (*l = IT + £i)$-„(z 2 = IT + £ 2 )) \ 1/X " _ 



\J{£\ +pi cos^i) 2 + (r + pi singly (£ 2 + P2 cos^) 2 + (r + p 2 sinfl 2 ) 2 

X 

(35) 



4pi P2 {£ 1 + Pl COS 9i ) (£2 + p2 COS 02 ) 

{£\ + £2 + Pi cos 0i + P2 cos 2 ) 2 + (pi sin X — p 2 sin 2 ) 2 



(£1 - (2 + Pi cos 0i - p 2 cos 2 ) 2 + (pi sin 6>i - p 2 sin 2 ) 2 ' 

The analytic continuation can be simplified for t, £%,£% S> e. For simplicity in the 
notation, but without loss of generality, we indicate as £1 the larger in absolute value 
of the two £ p and we assume £1 > 0. £2 can be either positive or negative. For t < \£ 2 \ 
we have a different result if the two points are on the same side or on different sides 
of the defect, namely: 

!log^P + 2ci 

S A {t<\£2\) = {l , e a _ (36) 



lQ g ;/ , / 2 2 + 2g i ^ 2>0 - 



c (ii - £ 2 y m^ 2 
e og (£1 + £ 2 y 

The first is just the sum of the entanglement entropies of two slits of length £\ and 
1^21 in the half-line starting at the origin. In fact, in the initial state the two parts 
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are unentangled and the total entropy is the sum of the two. Also the second result 
is just the ground state value for the slit [£2,£i] m the half-line [I]. 

The entanglement between A and B becomes sensitive to the joining of the 
system when the quasiparticles emitted from the origin arrive at £2- The behavior 
for \£ 2 \ < t < £1 does not depend on the relative sign of the £ p : 

SAm < t<i ^ = Q lQ S +t) ^ + 2ci ' (37) 

Finally for t > £\ we have almost the ground state value: 

S A {t>h) = l\J-±— e l + 2c' 1 , (38) 
6 a 

that is exactly the same as for case III. 

The resulting Sa for different values of £ p is plotted in Fig. [3] For \£ 2 \ < t < £\ 
there is again a large plateau whose actual extension and value depend on both £ p . 
The analytical value of the plateau is quite complicated and not really illuminating 
(it contains cubic roots, as can be realized solving the equation for the maximum). 
The most important features are that it decreases when the defect penetrates deeply 
in the slit (left part of Fig. \5§ and it stays almost constant when the slit moves far 
from the defect (right part of the figure). 

Finally one can easily read off from our formula the result for a defect exactly 
at the center of the slit: the entanglement entropy is independent of the time. This 
is different from what found numerically in a lattice model [13] and will be discussed 
(among the other things) in the next subsection. 



3.5. Comparison with numerical works 

As far as we are aware the entanglement entropy after a local quench has been discussed 
in only two papers [T21I2H]- There are a few qualitative differences between these results 
and the asymptotic ones of the CFT that are easily understood in terms of the general 
scenario we draw for quasiparticles emitted from the origin. These differences are 
attributable to slow quasiparticles, that exist as consequence of a non-linear dispersion 
relation Vk — dkEk, emitted from the joining point at t = 0. 

The first difference is that the lattice numerics present, on top of a smooth curve, 
very fast oscillations. These oscillations have been discussed in the context of global 
quenches [THIIH], and are due to the modes at the zone boundary \ka\ = tt that have 
Vk = 0. However, they are only corrections to the asymptotic result for t,£ a, since 
their amplitude remains constant while the CFT result diverges. Furthermore for any 
specific model they can be easily predicted in a quantitative way. 

In Ref. [13] different configurations of slits were considered for electron hopping 
on a chain (XX model in spin language). The sharp horizon effect is present even in 
the numerics, since there cannot be (by definition) quasiparticles faster than the CFT 
ones. We already discussed that the numerical results for a defect at the boundary 
of the slit is in quantitative agreement with the CFT in the region t < £. However 
for t > £ the numerics show an asymptotic value that appears to be exactly that of 
the ground state. Accordingly to our analysis this is possible only when the boundary 
entropy logg vanishes. 

The most relevant qualitative difference is that for t greater than the maximum 
length the asymptotic result is not constant but it is slowly decaying (like log(£)/£) 
smoothing the sharp transition from the plateau. This, for example, changes 
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completely the behavior in the case of central defect. We expect that this phenomenon 
can be ascribed to the slow quasiparticles that, after the quickest ones have arrived 
entangling the two parts, then disentangle them for very long times. This is the most 
plausible explanation, but we do not have an argument to put it on a more quantitative 
level. 



4. Correlation functions 



In this section we derive the time-dependence of correlation functions after the local 
quench. We find that the one-point function whose functional form is completely fixed 
by conformal invariance. The two-point functions instead depend on the particular 
(boundary) CFT. We will discuss all the details in the simplest case of a gaussian 
theory and then show how from general CFT arguments we can obtain part of the 
asymptotic behavior. Some of these results may have been previously derived in the 
context of quantum impurity problems from a different point of view (see e.g. |27j). 



One-point function 
The one-point function of a primary field in the half-plane Re w > is 

(*(">)) = roT^ix* » ( 39 ) 
|2Reiu| x * 

where a:$ is the scaling dimension of the field and Af is a non-universal amplitude 
that can be fixed in terms of the normalization of the two-point function of the same 
operator. We also fix the lattice spacing a to 1. Af is known for the simplest 
universality classes as the Gaussian theory and the Ising model 28j . 

With the mapping (jSJ) we can get the one-point correlation, that obviously 
assumes the same form as Eq. (|20[) . At the point r, after continuing to real time 
we have 

Aj(2r)- X * t<r, 

{<!>(,, n) { / £ y (40) 

Ab U* 2 - r 2 ) ) 

Thus for short times the correlation takes its initial value, until the effect of the joining 
arrives at time t = r when it decays for Or like t~ 2x<s > (note that this exponent is 
twice the boundary one). 



4-2. A two-point function: the gaussian model 

For the gaussian model the two-point function of a primary field <I> in the half-plane is 
given by Eq. (|55|) with — > (?~2,ri) and c„ — > {Af) 2 — 1 14J. As a consequence 
we can obtain its scaling as a byproduct of the result for the entanglement entropy 
for a general slit in Sec. 13.41 

For simplicity in the notation we assume here and in the following section n to 
be positive and to be larger than the absolute value of r 2 that can be either positive or 
negative. Thus we have that for t < 1 7*2 1 the two-point function keeps its initial value 
that has a different form depending on the relative signs of r p given by |4rir2| _a: ' 1 ' 
(different sides of the defect) and |4rir 2 (ri — r 2 ) 2 /( ?, i + r2) 2 \~ XiS ' (same side). 
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For t > T"i, it reaches the ground state value \r\ — r 2 \ 2x * ■ An interesting and 
non-trivial behavior is displayed for \r 2 \ < t < 7*1 when 



($(n,t)*(r2,t)) 



4-3. The general two-point function 



{ri+r 2 )(r 2 +t) 



(ri — r 2 ){ri — t) 4ri(i 2 — r. 



2\ 



(41) 



There are some features of the previous result that are expected to be valid in general, 
not only for a Gaussian theory as for example the the horizon effect and the final 
equilibrated value. The natural expectation is that these results are obtainable with 
CFT and in fact we will show that this is the case. 

The two-point function in the half-plane can be always written as |14j 

<*(«*)*(«*)> = f 1 ^ 1+ ?j^T?j, " F(r,) , (42) 

\ \Wi — W 2 \\W 2 — Wi\\Wi + Wi\\W 2 + w 2 \ J 

where 77 is the four-point ratio 

_ \Wi + Vh]\W2 +W 2 \ 

V 1 1 - 11 1 - 1 ' \^) 

\Wi + W 2 \\W 2 + I0i| 

and the function F(rf) depends explicitly on the considered (boundary) model. It 
is known for the simplest models as e.g. the Gaussian [F(rj) = 1) and the Ising 
universality class (see below). 

Thus the behavior of the two-point function depends mainly on the value on the 
ratio 77 that we now study. Using the analytic structure of the previous section we 
have for t < \r 2 \ a different behavior if the two points are on the same or on different 
sides of the defect: 

4rir 2 



(n + r 2 ) 2 



T2 > 0, 

r 2 < 0. 



For intermediate times we have 

v(w<t<r 1) = 1 2r : { vitv ^ 

{ri + r 2 )(ri + t) 

while for larger times 77 is 1 constantly. 

From the previous subsection we already know how the first part of Eq. (|42p 
transforms under the conformal mapping ([5]). Thus we need only to map F{rf) that, 
in general, is an unknown function. However we know its behavior in two special 
circumstances. Indeed when r\ ~ 1 the two points are deep in the bulk, meaning 
F(l) — 1. Instead for 77 <C 1, from the short-distance expansion, we have 

where Xb is the boundary scaling dimension of the leading boundary operator to which 
$ couples and Af is the bulk-boundary operator product expansion coefficient that 
equals the one introduced in Eq. (|3"5|) [see e.g. [2"5]]. 

From this we easily understand the behavior for t > r%. Since 77 — > 1 so F(X) = 1, 
the two-point function is just the one in the ground state 

($(n < t)$(N < i)> = (47) 

\Tl — 7-2 * 
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Although this result could have been expected, it is important to have recovered it 
only from CFT arguments. 

The behavior for |r 2 | < t < r\ is instead more complicated: combining the result 
of the previuos section with the value of r\ we have 



($(n,t)*(r2,t)) 



(n + r 2 )(r 2 + t) 



(n - r 2 ){n -t) 4n(t 2 -r 2 ) 



F 



2n(r 2 +t) 
(ri+r 2 )(n + t) 



(48) 



For example when $ is the order parameter in the Ising universality class we have [14 



V2 



(49) 



and x$ = 1/8. The sign ± depends on the boundary conditions. + corresponds to 
fixed boundary conditions and and — to free ones. This completely fixes the behavior 
of the two-point function in the intermediate regime. It would be very interesting to 
check this prediction. 

The time evolution for t < \r 2 \ depends on the sign of r 2 . For r 2 > we have 



($(n,t)$(r 2 ,t)) = 



1 {ri+r 2 f 



F 



4rir 2 

(n + r 2 f 



(50) 



_4rir 2 (n — r 2 ) 2 _ 

i.e. it keeps its inital value until t — r 2 , that is nothing but correct boundary value at 
zero time. 

More complicated is the behavior for r 2 < 0, in fact we get 



($(ri,t)$(ra,t)) 



1 



|4nr2| x<: 

$\2 



-F 



e 2 r 1 \r 2 \ 



{rf-t 2 )(r. 
e 2 r 1 \r 2 \ 



t 2 ) 



(rf 



t 2 )(r 2 2 



t 2 ) 



(51) 



|4rir 2 | x< " 

where in the last approximation we used that t <C t,r\,r 2 and the behavior of F(rj) 
for r] <C 1. In particular we note that this is time independent only when Xb = 0, 
which corresponds to the case of a non-zero 1-point function. Similar anomalous 
time-behavior was found when Xb > in Ref. [BJ. 



5. Decoupled finite interval 

A natural question arising is how the results we just derived change when we introduce 
more than one defect in the line. It is straightforward to have a path integral for the 
density matrix: we only need to have pairs of slits for — oo to —it and from it to 
-Moo everywhere there is a defect. However it becomes prohibitively difficult to treat 
this case analytically. In order to begin to understand the case when a finite interval 
interval is initially decoupled, we consider the case when it lies at the end of a half-line. 

So, let us consider a semi-infinite chain in which the A subsystem is the finite 
segment (—£, 0) and the B is the complement (0, oo) and with the initial defect at 
ru = 0. The space-time geometry describing this situation is like the one just 
considered, with a wall at — I + iy (y real) that represents the boundary condition. 
This is depicted in the left panel of Fig. |U 

In these circumstances the inverse conformal mapping between the z plane and 
the half-plane can be worked out using the Schwarz-Christoffel formula. After long 
algebra one obtains 

/ £ —{ w — l\ 

z(w)=i(-log(iw) + b— — — , (52) 

V 7T — IW +11 
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2e 





w 



Figure 4. Left: Space-time region for the density matrix of a semi-infinite system. 
Center: Mapping of the z-plane under the transformation Eq. J9]l- Right: Final 
map to the half-plane of the £ plane given by Eq. jSP when i 2> e. 



with the parameter b related to £ and e in a non-algebraic way. (A slit in the full line is 
closely related to this transformation, the last piece is replaced by (w 2 — l)/(w 2 + 1).) 
Unfortunately the mapping (I52p is not analytically invertible and its exact use is 
limited to numerical calculations that do not help us, since we need to perform an 
analytical continuation. We will develop in the remaining part of the section an 
approximate solution for < » 6, that however is not always justified: this limit is 
allowed only after the analytical continuation to real time. 

The limit £ 3> e before the analytical continuation simplifies the calculation 
because under the conformal transformation ((9|) the boundary at x = —£ is 
approximately a disk tangent to the imaginary axis at the origin and with radius 
R = e/4£ <C 1, as depicted in the central part of Fig. 0J This is simply checked 
applying the inverse transformation © to the disk £ = R(l + e 10 ) with 8 E [0, 2ir}. In 
fact we have 

i0)2( 1 + e -i0)_ ( 1 + 



C 2 -l R\l 



4R(1 



The half-plane minus the disk R(l 
transformation 



~ = ~4l?~ Z 4l? tan0/2+ ° (i?) ' (53) 



■ e ) is mapped in the half-plane Kew > by the 



w 



-iexp{2mR/0 . (54) 

Let us note that when e £, the parameter b is given by b = ne 2 / (8£) + 0(e 4 /£ 2 ). 

Combining the two transformations we can map the space-time region on the left 
of Fig. [3] into the half-plane (only in the limit £ 3> e) by means of 



-i exp 



me 
~2£ 



(55) 



with inverse 

£ . . . e 2 n 1 
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Note that the residue of the pole at w = — i is the same as the exact one in Eq. (|S"2"|) , 
Furthermore, calculating the the difference between the approximate and the exact 
solution, one easily checks that it is of the order of e 4 /£ 2 everywhere in the complex 
plane. 

We then have that w(ir) can be written as 



w(ir) 
Thus 



-te 



3 (vi- 



1 /e 2 Tve/2l) + i sin(V 1 - T 2 /e 2 7re/2f) 



($„(ir)) = c n 



t 2 sin(7rv? 



V2£) 



(57) 



(58) 



Note that for I — > oo it reduces to the previous result, as it should. Furthermore, it is 
clear that it is well defined only for I ^> e. 

Continuing to real time r — > it, and for f>e we have 

1 



($ n (ir)) = Cn 



TIC 



(59) 



_4^sin(7rf/2f)_ 

Clearly this cannot make sense when the argument of the power-law becomes negative 
(i.e. for t > 21), signaling that there is something wrong in the derivation. However, 
using the replica trick, for the entanglement entropy we obtain 



c ( M 
S A = - log[—tsin(Trt/2£) 

6 V 7T£ 



(60) 



One is tempted to assume that this result can be correct only for t < I and that 
for larger time it saturates as suggested by the quasiparticle interpretation. Only an 
exact calculation via the exact mapping (|52p can resolve these doubts. However exact 
solutions of integrable models or numerical density matrix renormalization group could 
already be able to eventually exclude Eq. (|60[) and to shed some light on the problem. 
We mainly mention this topic here to encourage further studies in this direction. 

Finally let us point out that transformations similar to Eq. (|52p appear when 
studying finite-size effects. 



6. Discussion 



We presented a detailed study of the unitary evolution that results after a local quench 
in a quantum one-dimensional system. All our findings have been obtained by means 
of CFT and so they are limited to the evolution of a gapless quantum model with a 
linear dispersion relation. We presented calculations for the entanglement entropies 
for different bipartitions of the system and for one- and two-point correlation functions 
in a general CFT. In particular we studied the case when two half-systems are joined 
at time t = 0. 

All our results are interpretable in terms of a scenario that we believe to be valid 
in general, not only for gapless systems (in analogy to the case of global quenches 
[51 Uni [H] ) . In fact, the initial state is expected to generate quasiparticle excitations 
at ro =0 that then propagate freely through the system and carry all the information 
about entanglement and correlations. In the case of a CFT all these excitations travel 
at the same speed v s (= 1 by normalization). However in general, as a consequence 
of a non-linear dispersion relation, a full spectrum of velocities is expected. 

Although the physical cut may seem a rather specific situation, we believe that 
our findings are quite general. For example, a general defect (e.g. a weakened bond) 
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is asymptotically equivalent to the results we obtained here as long as the defect is a 
relevant perturbation. This excludes the XX chain that is the only model studied so 
far [T3] , where the defect is believed to be marginal [25] , but should apply to XXZ spin- 
chains [30j . Furthermore we expect a qualitatively similar behavior when a system 
has been artificially prepared in a configuration that is only locally different from the 
ground state. For example there are already several studies where the local horizon 
effect is evident (see e.g. those considered in Refs. 31 ]). A full quantitative analysis of 
most of these settings should be possible properly adapting our CFT treatment (e.g. 
with the theory of boundary condition changing operators [24] ) . On the other hand, 
also the case of an initial state with one or more kinks displays a time-evolution with 
an horizon effect similar to the one considered here (see e.g. Refs. [32] where also 
finite temperature results are obtained in a rigorous manner). It would be interesting 
to understand whether these results can be recovered analytically continuing some 
CFT results. Work in this direction is in progress. 

We described here several different physical situations, but many problems are 
still left for future investigations. A first goal would be to check our predictions in 
exactly solvable models. In fact, as far as we are aware only the paper by Eisler and 
Peschel [13] discussed these topics in the context of the XX model and their results 
are fully compatible with ours. 

Another simple generalization of our results is to study gapped exactly solvable 
models, like the Ising model in a transverse magnetic field or others admitting a free- 
field representation. In these case it is also interesting to understand whether some 
predictions can be made on the basis of the generalized Gibbs ensemble 33 . 

A more difficult question to study analytically concerns the role played by 
quenched disorder. This is expected to change at a qualitative level the quasiparticle 
scenario as a consequence of Anderson localization [34]. In this case the time- 
dependent density matrix renormalization group [35] should be quite effective, since 
for clean systems we know that the entanglement entropy never increases dramatically 
as it does in the case of global quenches. 
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